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Radio supernovae (RSNe) are weak and rare events. Their typical maximum radio luminosities are of the order of only 
10"'' ergs"' Hz"'. There are, however, very few cases of relatively bright (and/or close) RSNe, from which the expansion of the 
shock and the radio light curves at several frequencies have been monitored covering several years. Applying the standard model 
of radio emission from supernovae, it is possible to relate the defining parameters of the modelled expansion curve to those of the 
modelled light curves in a simple algebraic way, by assuming an evolution law for the magnetic field and for the energy density of the 
population of synchrotron-emitting electrons. However, cooling mechanisms of the electrons may affect considerably this connection 
between light curves and expansion curve, and lead to wrong conclusions on the details of the electron acceleration and/or on the 
CSM radial density profile. In this paper, we study how electron cooling modifies the flux-density decay rate of RSNe for a set of 
plausible/realistic values of the magnetic field and for different expansion regimes. We use these results to estimate the magnetic 
fields of different RSNe observed to date and compare them to those obtained by assuming energy equipartition between particles 
and magnetic fields. For some of the best monitored RSNe, for which deceleration measurements, optically thin spectral index, and 
power-law time decay have been observed (SN 1979C, SN 1986J, SN 1993J, and SN2008iz), we find self-consistent solutions for the 
index of the power-law circumstellar density profile (s = 2 for all cases), the index of the power-law relativistic electron population 
(rather steep values, p = 2.3 - 3.0) and the initial magnetic field (ranging from ~ 20 to > 100 G). 
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1. Introduction 

Radio supernovae (RSNe), which are the radio counterparts of 
core-collapse supernovae (SNe), are weak and rare events. Only 
about 10 - 20% of the observed SNe are detected in radio (e.g., 
Weiler et al. 2002). Moreover, their typical maximum radio lu- 
minosities are of the order of 10^^ erg s"' Hz"' (flux densities of 
the order of 1 mJy for extragalactic distances, close to the sensi- 
tivity limits of present detectors). There are, however, very few 
cases of relatively bright RSNe, from which the expansion curve 
of the shock, using Very Long Baseline Interferometry (VLBI) 
observations, and radio light curves at several frequencies were 
obtained covering, in some cases, several years, e.g.: SN 1979C, 
SN1986J, SN1993J, and SN2008iz. Although there ai-e only 
a handful of objects, their detailed study allowed to check and 
refine the current theoretical models of radio emission in super- 
novae. This small number of well-observed RSNe may also dra- 
matically increase in the near future, thanks to the forthcoming 
ultra-sensitive interferometers with a high spatial resolution, like 
the Square Kilometre Array (SKA). 

Using the standard model of radio emission from supernovae 
(Chevalier ll982all982bl) . it is possible to relate the defining pa- 
rameters of the modelled expansion curve to those of the mod- 
elled light curves in a simple algebraic way, by assuming an 
evolution law for the magnetic field (and for the density of the 
population of synchrotron-emitting electrons) and a radial den- 
sity profile for the circumstellar medium, CSM, (see, e.g. Weiler 
et al. 120021) . The decay in the radio-light curves according to 
this model is related to the time decay in the magnetic field and 
the radial decay of CSM density. However, the continuous en- 



ergy loss by the relativistic electrons (electron cooling), mainly 
due to synchrotron radiation (i.e., radiative cooling), but also 
to adiabatic expansion and inverse-Compton scattering, are not 
considered in the derivation of this relationship between light 
curves and expansion curve. Electron cooling may affect consid- 
erably the shape of the light curves for a given expansion curve. 
For instance, Marti- Vidal et al. ( 1201 Ibl) succesfuUy modelled the 
exponential-like decay of the SN 1993 J radio light curves at late 
epochs, reported in Weiler et al. (120071 ). using only radiative- 
cooling effects, and assuming that the density of the CSM was 
negligibly small at large distances to the progenitor stafl. In any 
case, it seems clear that if electron cooling is not considered in 
the modelling of the radio light curves of a supernova, it could 
result into wrong estimates of the model parameters. In this pa- 
per, we study how electron cooling modifies the flux-density de- 
cay rate of RSNe for several values of the magnetic field and for 
different expansion regimes. These results can be used to esti- 
mate the magnetic fields of observed RSNe. 

In the next section, we outline the standard model of ra- 
dio emission from supernovae. In Sect. [3] we study the effect 
of electron cooling in the population of emitting electrons and 
in the flux-density decay rate. In Sect. |4] we present the results 
of several simulations of the expansion and radio light curves 
of RSNe. In Sect. |5] we explain how these results can be used 
in real cases to estimate physical quantities in RSNe and esti- 



' An additional effect due to the escaping of the electrons from the 
emitting region might also be necessary to model the light curves of 
SN 1993 J, were the density of the CSM not negligible at those large 
distances to the progenitor 
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mate the magnetic fields for several observed RSNe, comparing 
these estimates to those obtained by assuming particle-field en- 
ergy equipartition. In Sect. |6] we summarize our conclusions. 



2. Connection between expansion and radio light 
curves in RSNe 

In the standard model of emission from supernovae (Chevalier 
11982 a, 1982b), the spherically-symmetric expanding shock is 
described as a contact discontinuity plus two shocks, one mov- 
ing backwards (from a Lagrangian point of view) and the other 
moving forward, shocking the CSM. A fraction of shocked CSM 
electrons is accelerated to relativistic energies, possibly due to 
statistical Fermi processes, and produce synchrotron emission at 
radio wavelengths as they interact with high magnetic fields in 
the shocked CSM region. 

The distance, r, from the contact discontinuity to the center 
of the expansion evolves as a power-law of time (r oc f'") with an 
expansion index, m, that depends on the radial density profiles of 
CSM (pcsM r'") and ejecta (pej oc r"") in the form (Chevalier 
ll982al[T9^ 



n-3 



(1) 



This solution of the shock expansion holds for n > 5 and 
s < 3. The structure of the shock (contact discontinuity plus 
backward and forward shocks) expands in a self-similar way. 
Therefore, the expansion of the forward and backward shocks 
also follows the law oc f'". 

On the other hand, the distribution of relativistic electrons in 
energy space follows a power law (A^ oc E^'') and the energy- 
density of the magnetic field is assumed to be proportional to 
the energy-density of the shock (i.e., oc n V^, where is 
the average magnetic field squared, n oc r^^ is the particle num- 
ber density, and V oc /"'-i)/"' jg ^jjg shock expansion velocity). 
Hence, 



Boct' 



,m(2-s)/2-l 



(2) 



We must notice a limitation in the standard model at this 
point. For a standard CSM particle density of 10*^ cm"^ at a dis- 
tance of 10'^ cm from the explosion center, and an expansion 
velocity of 20 000 km s"', a magnetic field of 50-60 G translates 
into a similar energy density for the expanding shock and the 
magnetic field. Such a large magnetic-field energy density may 
affect the hydrodynamics of the shoclfl This effect is neglected 
in the model (which, indeed, assumes that the magnetic-field en- 
ergy density is a small fraction of that of the shock). Hence, for 
cases of very large magnetic fields reported in Sect. 5, high CSM 
particle densities and/or large expansion velocities might be ac- 
cordingly considered, to make the magnetic-field estimates con- 
sistent in the frame of the standard model. 

The fraction of accelerated particles by the shock, or injec- 
tion efficiency of the shock, is also assumed to be proportional 
to the shock energy density. Under all these assumptions, and 
considering that the intensity of synchrotron radiation is (e.g. 
PachoIczvk [T970] l 



IccNB' 



,(l+p)/2 



(3) 



^ Detailed magneto-hydrodynamic simulations would be necessary 
to study the real impact of large magnetic fields in the evolution of the 
expanding shock 



it is possible to derive the intensity, /, in the optically-thin regime 
if we neglect electron cooling. Since, in that case, M(E^ oc 
E^P nr^V df, it can be shown that / oc v " f, with 



P ■ 



1 



and 



P^-{m{2{n+p)-{5+p)s)-2p-lO). 



(4) 



(5) 



This equation brings a direct relation between the decay in- 
dex of the radio light curves in their optically-thin regime, (3, on 
one hand, and the supernova expansion index, m, the energy in- 
dex of the injected relativistic electrons, p, and the index of the 
CSM radial density profile, s, on the other hand. For the case of a 
constant pre-supernova mass-loss wind (i.e., s -2) this equation 
reduces to/3 = (6m - p-5)/2 (e.g., Weiler et al. l2002b . 

3. Radiative and adiabatic cooling of the relativistic 
electrons 

The supernova shock is continuously accelerating electrons from 
the shocked CSM. These electrons are distributed as oc £ ''. 
However, the electrons already shocked that are emitting syn- 
chrotron radiation loose energy and, therefore, shift towards 
lower energies in the electron-energy distribution. Since the 
number of electrons is conserved, we can make use of the conti- 
nuity equation in energy space, i.e.. 



N = Ve(NE) + S(E, t) - HE, t). 



(6) 



where S (E, t) is the source function (the new electrons contin- 
uously accelerated by the shock) and L{E, t) accounts for the 
escaping of electrons from the emitting region. We will assume 
that L{E, = (in Marti- Vidal et aI. I20IIbl we use L(£, f) oc 
to model the SN 1993 J radio data, although the effects of this 
term are very small compared to S (E, t) until very late epochs, 
when a large drop in the CSM density profile takes place). It can 
be shown (see Appendix|A]l that the source function is S it) E^'', 
where 



S{t)=NoF,., 



1 



m(5-j)-3 



(7) 



where A^o is the number density of shocked CSM electrons at a 
reference epoch (fo), Frei is the fraction of accelerated electrons 
(of the order of 10"^ for SN 1993 J), and £■„, is the minimum 
energy of the relativistic electrons (we set E„ = c^, although 
this value is not relevant in the optically-thin regime of the light 
curves). 

The term E takes into account the energy loss (or gain) of 
the electrons. The energy loss can be either radiative, adiabatic, 
and/or due to free-free interactions with atoms or ions in the 
CSM. The energy gain can be due to self-absorption of the syn- 
chrotron radiation or to inverse-Compton scattering, although 
these effects are negligible in the optically-thin part of the light 
curve (and also for large magnetic fields), which is that of our 
interest here. In the case of radiative losses, we have 



E,^-C2BiE'\- 



(8) 



where C2 = 2.37 x 10 ^ in cgs units (see Pacholczvk l 1 9701 ), and 
is the magnetic field at a reference epoch (fo) averaged in 
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the orthogonal planes to the electron trajectories. For a random 
distribution of magnetic-field lines and electron trajectories, 
is equal to V2/3 times the total averaged magnetic field, Bq, at 
the reference epoch. In the case of adiabatic losses, we have 



Idr E 

E — m—. 

r dt t 



(9) 



Therefore, if radiative cooling and adiabatic expansion are 
the dominant processes of energy loss by the electrons, we have 



E = -C2B] E^ I - 



m{2-s)-2 



E 

— m—, 
t 



(10) 



In Eq. (To] we have neglected the term due to free-free in- 
teractions of the electrons with the surrounding CSM atoms and 
ions {E cc r"' E), since this term is much smaller than the ra- 
diative and adiabatic terms in the optically-thin regime of the 
light curves. In AppendixiBl we analyze under which conditions 
might the free-free term not be negligible compared to the radia- 
tive and adiabatic terms. 

Equation |6l together with Eqs. [T] and [TOl is a typical 
difussion-like partial differential equation that can be numer- 
ically integrated using, for instance, a semi-implicit approach 
(e.g. Marti- Vidal et al. 201 lb). However, since synchrotron self- 
absorption, inverse Compton, and free-free interactions are ne- 
glected (i.e., only the radiative and adiabatic terms in E are con- 
sidered), it is also possible to find an integral form for the so- 
lution of this simplified version of Eq. |6] We show this solution 
in Appendix ICl From the numerical solution of N(E, t), we can 
estimate the flux-density decay rate of the light curves, since the 
intensity is 



/oc f' 



m{2-s)i2-\ 



I 



NF{x)dE, 



where the power-law of time is related to the decay of the mag- 
netic field (see Eq. |2]i, x is the ratio between the observing fre- 
quency and the critical frequency at energy E, and Fix) is 



Fix) 



Ksjiiz) dz. 



being K^/^iz) a Bessel function of the second kind (e.g., 
Pacholczvk 1 1 9701 ) . Then, from the time evolution of /, we can 
estimate /3 for different combinations of m, p, s, Eq, No, and 
and compare the results to Eq. |5]in order to check the effect of 
electron cooling in the light curves. 

4. Effect of magnetic fields in the radio light curves 

We show in Figs.[T]and|2]the jS obtained from our simulations as 
a function of m and Bq (the magnetic field at the reference epoch 
f() - 5 days) for 6 values of p (2.0, 2.2, and 2.4, in Fig.[Tl 2.6, 2.8, 
and 3.0, in Fig.|2]l and for 3 values of s (1.6, 2.0, and 2.4). We 
have computed at 5 GHz between 300 and 1000 days after the 
shock breakout. Different selections of frequencies and/or age 
ranges result into deviations in of a few % at most. 

Since /3 is computed in the optically-thin part of the radio 
light curves, A^o is not really important in the simulations (chang- 
ing this value would affect the opacity in the early supernova 
evolution). In our case, the important quantity would be A^o F,-ei, 
which accounts for the number of relativistic electrons. Indeed, 
A^(£', t) only depends on Bq regardless of a constant scaling fac- 
tor defined by A^o Fyei- Therefore, the value of A^o Frei does not 



really affect the estimates of [5. To ensure that this statement is 
correct, we checked that the values of fi derived from our simu- 
lations are only sensitive to changes in Bq, m, p, and s. 

If radiative cooling is negligible (i.e., for small values of 
Bq), the j6 computed from our simulations approaches the values 
computed from Eq. |5]for all combinations of m, p, and s. This 
is an expected result, since the adiabatic losses alone (which are 
oc E) do not affect the power law of the electron distribution 
(Pacholczyk l 19701 1. However, as the magnetic field increases, /3 
decreases in absolute value (i.e., the light curves become flatter). 
This result is in principle non-intuitive, since one would expect 
the light curves to be steeper as the radiative cooling (i.e., the en- 
ergy loss of the electrons) is more important. The light curves in 
the optically-thin stage are flatter for larger Bq, because the mag- 
netic field decreases as the supernova expands (see Eq. |2]i and, 
therefore, cooling effects (which are smaller for smaller mag- 
netic fields) are less important as time goes by. Thus, A^ assymp- 
totically approaches the value without cooling as the supernova 
expands. As a consequence, N/N (which affects the value of /3) 
takes a larger value if we consider radiative cooling. In Appendix 
iDl we show the details of this discussion mathematically. 

The largest deviations of with respect to the cooling-free 
value (i.e., that of Eq.|5]l correspond, in all cases, to the smallest 
decelerations of the shock (i.e., values of m close to 1) and/or to 
the steepest CSM radial density profiles (i.e., larger values of s). 

In the case i = 2, we can approximate the /3 shown in Figs.[T| 
and |2] with the phenomenological equation 



y8(B()) 



Fi/il-m) 



/3iBo = Q) \Bo + Fi/il-m) 



F2/((l-m)(2p-l)) 



(11) 



where jS(Bo) corresponds to a magnetic field Bq and y6(B() = 0) 
is that given in Eq.|5](i.e., with no radiative cooling considered). 
The parameters Fi and F2 take the values 7.725 G and 0.184, 
respectively. The maximum deviation between the /? computed 
from Eq.fTTIand those shown in Figs.[T|and|2](for i = 2) is only 
3.5%. 



4.1. Changes in the spectral index 

If electron cooling is not considered, there is a direct relation- 
ship between p and the spectral index a (see, Eq.|4|: p - I + 2a. 
However, when electron cooling is taken into account, there is 
a flux of electrons towards smaller E, which increases the value 
of a. This effect is more important as we increase the observ- 
ing frequency. We must notice, however, that new electrons are 
continuously being injected in the emitting region, and their en- 
ergy distribution is assumed to be always oc E^p, so this fraction 
of electrons is not affected by cooling. Therefore, the effect of 
cooling in the spectral-index steepening is somewhat mimicked 
by the new electrons entering the shocked CSM. The integra- 
tion of Eq. |6] takes into account this trade-off between electron 
cooling and the source function. In Figs. [3] and H] we show the 
simulated spectral indices, averaged between 300 and 1000 days 
after shock breakout and centered at 5 GHz. We show a for the 
same values of p and s used in Figs.[T]and|2l 

We notice that radiative cooling is more important at higher 
energies, so the (effective) spectral index should slightly increase 
with the observing frequency, for instance, the difference be- 
tween the spectral indices centered at 5 GHz (which are higher) 
and those at 1.7 GHz (which are lower) is typically 2-3% for 
magnetic fields of 10 G and 5-6% for magnetic fields of 100 G. 

The values of a obtained from our simulations tend to the 
expected values without cooling (i.e. a - ip - l)/2) when the 



4 I. Martf-Vidal et al.: Electron cooling in radio supernovae 



s = 1.6 ; p = 2.0 s = 2.0 ; p = 2.0 s = 2.4 ; p = 2.0 




Fig. 1. y6 as a function of m and Bo (reference epoch fo = 5 days) for a selection of values for s and p. For Bo ~ 0, we roughly recover 
the yS given in Eq. |5] As we increase Bq, /3 approaches 0. This effect is more pronounced for larger values of s (i.e., for steeper CSM 
radial density profiles). 
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s = 1.6 ; p = 2.6 s = 2.0 ; p = 2.6 s = 2.4 ; p = 2.6 




-1.35 -1.20 -1.05 -0.90 -0.75 -0.60 -0.45 -0.30 -0.15 -1.90 -1.70 -1.50 -1.30 -1.10 -0.90 -2.48 -2.32 -2.16 -2.00 -1.84 -1.68 



/3 /3 /3 



Fig. 2. Same as Fig.lT] but for a different set of values for p. 



magnetic field approaches 0, also as expected. An increase in Bo steepens the spectrum of the radiation (i.e., a increases), for 

combinations of s, m, and p, because Z^rad oc E^. 
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s = 1.6 ; p = 2.0 s = 2.0 ; p = 2.0 s = 2.4 ; p = 2.0 




Fig. 3. Spectral index a, centered at 5 GHz, as a function of m and Bq (reference epoch to = 5 days) for a selection of values for s 
and p. For Bq ~ Owe roughly obtain the canonical value a = (p - 1)/2. As we increase Bq, the spectra become steeper (a increases). 
This effect is more pronounced for lower values of s (i.e., for flatter CSM radial density profiles). 




s = 1.6 ; p = 2.8 



s = 2.0 ; p = 2.8 




20 40 60 80 100 





0.90 0.95 1.00 1.05 1.10 1.15 1.20 1.25 1.30 1.35 0.88 0.96 1.04 1.12 1.20 1.28 0.90 0.95 1.00 1.05 1.10 1.15 



s = 1.6 ; p = 3.0 



s = 2.0 ; p = 3.0 



s = 2.4 ; p = 3.0 




1.00 1.08 1.16 1.24 1.32 1.40 1.00 1.04 1.08 1.12 1.16 1.20 1.24 1.28 1.32 1.36 1.00 1.05 1.10 1.15 1.20 1.25 



Fig. 4. Same as in Fig. [3] but for a different set of values for p. 



In the case s = 2, we can also approximate the a shown in where a(Bo) corresponds to a magnetic field Bo and a(Bo = 0) is 
Figs. [3] and ID with a phenomenological equation (p - l)/2 (i.e., no radiative cooling considered). The parameter 

+ 11 , (12) 



a(Bo) 



aiBo^O) \Fi(2m-l)(2p + 5) 



8 



I. Martf-Vidal et al.: Electron cooling in radio supernovae 



Fi takes the value 3.04 G, and the maximum deviation between 
the a computed from Eq. [12] and those from the simulations is 
also 3.5%, for all the analyzed values of m, p, and Bo- 
Figures [TJIJIEI and|4](and eventually Eqs.[TT]and[T2]i can be 
used to estimate the magnetic field in a supernova by using the 
a, m, and /3 inferred from the observations (provided light curves 
and the expansion curve of the supernova have been observed). 
In the next section, we will estimate magnetic fields in a number 
of radio supernovae, for which expansion curve and radio light 
curves are available. 

4.2. The special case p-2 

We notice that for the special case p-2, the effect of radiative 
cooling in the electron energy distribution should be negligible 
for all E, since NE^ would not depend on E and its contribu- 
tion to the energy gradient of would therefore be null (see Eq. 
|6]). In principle, one would expect the population of electrons to 
evolve as if there were only adiabatic cooling, so neither fi nor a 
should depend on the magnetic field. However, when p-2, the 
evolution of N{E, t) is not only determined by the source func- 
tion, S (E, t), and the adiabatic term, but also by tp (see Appendix 
ICl Eq. IC.6I 1. which is the time at which all the electrons with en- 
ergies larger than E have energies below E at time t. The time 
tf is larger than fo for high energies and/or large t, and depends 
on the magnetic field. Thus, even for p-2, the light curves and 
spectra will be somewhat modified by radiative electron cooling 
at high frequencies and late epochs (those frequencies and su- 
pernova ages depend, of course, on the strength of the magnetic 
field and the deceleration of the shock), as it is shown in Figs.[T] 
and |3] (upper rows) . 



5. Estimate of magnetic fields in observed RSNe 

If a radio supernova is strong enough to be monitored with 
VLBI, it is possible to estimate m from the expansion curve and 
P and a from the light curvefl If cooling is not considered, from 
a it is possible to derive p and, using Eq. |5] it is possible to 
derive s. Additionally, assuming a constant temperature of the 
CSM electrons, the opacity due to free-free processes decreases 
as f, being 6 - m{\- 2s) (e.g., Weiler et al. 2002). Therefore, if 
the light curves are well sampled in the optically-thick regime, 
another condition can be imposed to the parameters if we as- 
sume dominance of free-free absorption (FFA) over synchrotron 
self absorption (SSA). Self-consistency between all the parame- 
ters can then be checked. 

However, it is not clear how much FFA dominates the light 
curves of usual RSNe. For instance, SSA has shown to be, by 
far, the dominant absorption mechanism in all the evolution of 
the SN 1993J light curves (Fransson & Bjornsson 1998; Perez- 
Torres et al. 2001 ; Marti- Vidal et al.^OHb). Moreover, different 
forms of electrons cooling, as we have shown in the previous 
section, affect the values of yS and a for a set of m, s, and p, 
depending on the strength of the amplified magnetic field. In this 
section, we will infer the values of magnetic fields of several 
RSNe, based on the the values of a, /3 and m estimated from the 



' We assume that all these quantities are determined in the optically- 
thin regime, which corresponds to a positive a (i.e., a steep spectrum) 
and a decreasing (or non-increasing) flux density, with the exception of 
very special cases (s « 2 together with m ~ I, see Figs. 1 and 2). In 
all the observational cases studied in this paper, the conditions for an 
optically-thin regime hold for all the quantities used in our analysis. 



observations. An a priori value for s and/or p must be however 
assumed to estimate Bo using this approach. 

5.1. SN1979C 

Weiler et al. (|1991| l reported more than 10 years of observations 
of the SN 1979C radio Ught curves at 15, 5, and 1.4 GHz. These 
authors fit a = 0.74+J5 JJ5 ^nd /3 = -0.78+° °2 regard of the 
expansion curve, different results have been reported by differ- 
ent authors. Marcaide et al. (I2002I I reported a strong deceler- 
ation in SN 1979C and, from a more complete VLBI dataset, 
Bartel & Bietenholz (120031 1 reported an almost free expansion 
(i.e., m ~ I) for 22 years. More recently, Marcaide et al. ( I2009bl ) 
re-analyzed their VLBI data and complemented them with new 
1.6 GHz observations and the data from Bartel et al. These au- 
thors arrived to the conclusion that the expansion of SN 1979C 
was indeed almost free (m - 0.91 ± 0.09) for over 25 years. 

The fitted a is very close to the value corresponding to p = 
2.5 without radiative cooling. Therefore, we conclude that either 
the magnetic field is very small (and hence a ~ (p - l)/2), or 
p is lower than 2.5. Assuming p = 2.2 (or p = 2.4) and s - 2, 
we estimate from Fig.[T]a magnetic field of ~ 20G (or ~ 40 G) 
at day 5. There are no solutions neither for 5 = 1.6 nor s = 2.4. 
Now, from Fig. [3] the observed a and m imply Bo ~ 20 - 30 G 
(for p = 2.2) and Bq ~ 5 - lOG (for p - 2.4). Therefore, based 
on the radiative-cooling assumption, the magnetic field at day 5 
should be between 20 and 30G if /? ~ 2.2. Indeed, from Eqs.[T2l 
and[TT]we find self-consistent estimates of a and /3 for p = 2.3 
and Bo = 30 G. 

How do these estimates ob Bo compare to the equipartition 
magnetic field? In the case of energy equipartition between par- 
ticles and fields, it is possible to estimate the magnetic field in 
the radiating region provided the size and the total luminosity of 
the source are known. The expression used for this estimate is 
taken from Pacholczyk (I1970I I: 

Be, = (4.5ci2(l + k)l4>fl'R-''^'Lf, (13) 

where c\2 depends on the spectral index, a, and on the frequency 
range in the spectrum integration. is the filling factor of the 
emitting region to a sphere (0.66 for a shell-like structure of 30% 
fractional width), R is the source radius, Lr is the integrated ra- 
dio luminosity, and k is the ratio between the heavy particle en- 
ergy density to the electron energy density. We do not know the 
details of the particle acceleration, and the efficiency of acceler- 
ation could depend on the particle mass. Hence, k can vary from 
1 (case of a much larger acceleration efficiency of the electrons 
compared to the protons) to njp/mf ~ 2 x lO"* (case of a similar 
acceleration efficiency for electrons and protons). 

Using the spectral index and flux densities by Weiler et al. 
(119911) . the expansion curve by Marcaide et al. (2009bll, and the 
distance to the host galaxy (M 100) by Ferrarese et al. d 19961 1 of 
16.1 Mpc, we estimate L« = 1.6 x lO"'-' ergs"' at day 5 after ex- 
plosion. Therefore, the equipartition magnetic field could range 
between 10 and 85 G (for k-\ and k - 2000, respectively). Our 
estimated Bo, assuming p - 2.3 and s - 2, corresponds to low- 
to-intermediate values of k, i.e. low-to-intermediate efficiency in 
the acceleration of ions. 

We must notice that the cooling-free prediction of fi (Eq. 
|5]l for s - 2 and p - 2.5 is consistent with the observed 
one if m - 0.99 (i.e. for an essentially non-decelerated expan- 
sion), which is indeed compatible with the value of m reported 
in Marcaide et al. (.2009bi) at a Icr level. This latter possibility 
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Table 1. Model parameters for several RSNe. Observed refers to 
those obtained from the fitted expansion and radio light curves; 
Assumed and Derived refer to those obtained from comparison 
with the results shown in Figs. [T][2l[3] and|4] Cases with two 
possible solutions are given in two rows, one row for each solu- 
tion. 



Supernova 




Observed 




Assumed 


Derived 




m 


a 






P 


So(G) 


SN 1979C 


0.91 


0.74 


-0.78 


2.0 


2.5 


~ 










2.0 


2.3 


20-30 


SN 1986J 


0.69 


0.7 


-1.18 


1.7 


2.4 


~ 










2.0 


2.0 


30-50 


SN 1993J° 


0.87 


0.98 


-0.78 


2.0 


2.5 


60-80 


SN 2008iz 


0.89^ 


1.08 


-1.43 


2 


~ 3 


~ 










2.4 


2.6 


~ 100 



" Fit to data between days 300 and 1000 after explosion (see text). 
* Derived from VLBI observations (Brunthaler et al. 120 lib 



would imply a very small magnetic-field energy density, com- 
pared to the energy density of the particles. 

In Table l5.11 (row 1) we summarize the values of s, p, and Bo 
discussed for this supernova. 

5.2. SN1986J 

There are a number of peculiarities in the radio light curves of 
SN 1986 J compared to those of other RSNe. The spectral in- 
dex cannot be well fitted to a single value from 1 .4 to 23 GHz 
(Weiler, Panagia & Sramek ll9901 l. In the optically-thin part of 
the light curves, a = 0.7 + 0.1 between 5 and 15 GHz, but 
a = 0.2 ± 0.2 between 15 and 23 GHz. Additionally, Bietenholz, 
Bartel, & Rupen (2004) reported the discovery of a compact 
source in the shell center of SN 1986J with an inverted spec- 
trum, and interpreted this source as due to accretion onto a black 
hole or to a young pulsar nebula. 

The best-fit parameters reported in Weiler, Panagia & 
Sramek (1990) ai-e a = 0.67+°°^ and /3 = -l.lS+JSJj^, based 
on observations up to year 1989. Bietenholz, Bartel, & Rupen 
(I2OO2I 1 reported a much lower /3 for later epochs that slightly de- 
pends on the observing frequency (ranging from -2.7 at 8.4 GHz 
to -3.5 at 23 GHz). In this work we will use the a and /3 obtained 
from the early epochs (i.e., those up to year 1989) and between 
5 and 15 GHz. 

In regard of the expansion curve, Bietenholz et al. ( 1201 Oi l 
reported m = 0.69 ± 0.03, a value much lower than those of the 
other RSNe observed with VLBI (-0.8-0.9). 

Now, from the extrapolated size at day 5, a distance to the 
host galaxy (NGC 891) of 8.4 + 0.5 Mpc (Tonry et al.'200TT), and 
using a - 0.7, we obtain an equipartition magnetic field between 
14 and 100 G (fc = 1 and 2000, respectively) using Eq.[T3]for day 
5 after explosion. 

How do these estimates of Bq compare to those that can be 
obtained with our approach? A spectral index of 0.7 can only be 
obtained with p - 2.4 or lower Trying with the lowest value, 
p = 2, we find from Eqs. [TT] and [12] (or Figs. [1] and |3]l self- 
consistent values of a and /3 with Bq - 30- 50 G. Using p = 2.2 
and i = 2, we estimate from Fig.[T]a magnetic field of ~60G 
at 5 days after explosion. Now, from Fig. |3] and assuming the 
same values for s and p, a magnetic field of ~10G is estimated. 
Both estimates are incompatible. For s - 2.4 or s - 1.6 we can 
neither obtain a coherent estimate of the magnetic field; using 



now p - 2.4, the observed a requires, of course, Bo ~ and 
the observed /3 can only be explained with our simulations if 
s ~ 1.6. Therefore, a compatibility between Figs. [T] and |3] is 
found for small values of the magnetic field and a rather flat 
CSM radial density profile (s ~ 1.6). Indeed, the cooling-free 
prediction of p given by Eq. |5] (which is similar to that one with 
cooling considered if Bq is very small) is equal to the observed 
one for p = 2.4 and s ~ 1.7. Hence, we conclude that either 
/9 = 2, 5 = 2, and Bo = 30 - 50G, or 5 < 2, p ~ 2.4 and 
Bo ~ G, can explain the radio data for this supernova. In table 
15.11 (row 2) we summarize the values of s, p, and Bo discussed 
for this supernova. 

5.3. SN1993J 

This is the radio supernova with best-observed light curves and 
expansion curve (see Perez-Torres, Alberdi & Marcaide 200JJ 
Bartel et al. f 2002l Marcaide et al. |2009at Weiler et al. .20071 
Martf-Vidal et aI. I201 lall201 Ibt and references therein). 

Fitting their observed light curves (taken until ~4900 days 
after explosion) Weiler et al. (120071) obtained a = 0.81, S = 
-1.88, and p = -0.73. Therefore, without considering electron 
cooling, from the fitted a we obtain p = 2.6 and, from the expan- 
sion index reported in Martf-Vidal et al. (2011b) at late epochs 
(m - 0.87), we obtain s - 1.6. Applying now Eq.|5] we derive 
j6 = -0.44, which is inconsistent with the value fitted to the light 
curves (/3 = -0.73). 

However, if we decrease m down to 0.82, we can obtain a 
self-consistent solution for /3, using Eq. |5] This seems to be a 
strong evidence of a CSM radial density profile with an index 
s < 2. Also, Mioduszewski, Dwarkadas & Ball (2001 ) simulated 
radio images and the radio light curves of SN 1993 J without tak- 
ing radiative cooling into account, and claimed that s ~ 1 .7 pro- 
vides the best fit to the data. 

However, the evidence of < 2 coming from Eq. |5j and 
from the fit of the optically-thick part of the radio lightcurves, 
holds as long as the temperature of the thermal CSM electrons 
is taken constant throughout the whole extent of the CSM (to be 
able to use 6 - m(l - 2s)), which is not likely to apply in the 
ca se of SN 1993J (Fransson & Bjornsson 1998, Maitf-Vidal et 
al. l201 Ibl) . Additionally, more recent analyses of the X-ray data 
from SN 1993 J also discard the models with s < 2 (Nymark, 
Chandra & Fransson 2009; Chandra et al. 2009). 

From their simultaneous analysis of the complete light 
curves and expansion curve of SN 1993 J, Martf-Vidal et al. 
dMIIB reported Bo = 65. 1 + 1.6 G and/? = 2.59 ±0.01 for 5 = 2. 
It was also noted by these authors that using values of s < 2 re- 
sulted in poor fits to the data. From a much time-limited set of 
flux-density measurements, Fransson & Bjornsson ( 1998 ) fitted 
a similar magnetic field for day 5 after explosion (Bq ~ 68 G) 
using also s = 2, although they fitted a different energy index for 
the electron distribution (p = 2.1). 

Which magnetic field do we obtain for SN 1993 J with our 
approach? Opacity effects in the supernova ejecta may affect 
the spectral index and /3 at different frequencies and for dif- 
ferent times (Marcaide et al. '2009a', Martf-Vidal et al. 1201 Ibl l. 
Therefore, in our approach we must use the values of yS and 
a fitted to the subset of data where such ejecta-opacity effects 
are minimum or non-existent, and not those fitted to the whole 
dataset. Using the 5 GHz and 8.4 GHz data of Weiler et al. 
(120071) from day 300 to day 1000 after explosion, we obtain 
/3 = -0.78 + 0.05 at 5 GHz and /3 = -0.79 + 0.08 at 8.4 GHz. The 
average spectral index between 8.4 and 5 GHz at these epochs is 
» = 0.98 ±0.19. 
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This spectral index implies p - 3.0 or lower. For any value 
of p, neither ^ = 1.6 nor s - 2A yield self-consistent estimates 
of Bo using our approach. This is an additional evidence of a 
CSM with 5 = 2 for SN 1993J. Assuming now that s - 2, we 
estimate from Eqs.[TT]and[T2]that p ~ 2.4 and Bq = 60 - 80 G. 

In table 15. II (row 3) we summarize the values of s, p, and Bq 
discussed for this supernova. 

The range of values of Bq estimated this way is in agreement 
with the estimates reported in Fransson & Bjornsson d 19981 1 and 
Marti- Vidal et al. (201 IbV Nevertheless, here we have used a 
subset of the observed light curves, to avoid the undesired con- 
tribution of ejecta-opacity effects in our rough radiative-cooling 
model. 

Fransson & Bjornsson d 19981 1 and Marti- Vidal et al. (1201 Ibl l 
also discussed on the particle-field energy equipartition, based 
on their fitted magnetic fields. In both papers, it is concluded 
that, to obtain energy equipartition, an acceleration efficiency of 
the ions similar to that of the electrons (i.e., A; » 1 in Eq. [TlT l 
should take place in the shock. 

5.4. SN2008iz 

Marchili et al. ( 120101) reported a 5 GHz light curve for this super- 
nova, taken with the Urumqi telescope. Brunthaler et al. ( 1201 Oi l 
reported VLBI observations from which the explosion date and 
the expansion velocity could be estimated. Marchili et al. (1201 Oi l 
estimated an equipartition magnetic field between 0.3 G and 
2.1 G (for k - 1 and k = 2000, respectively) at day 63 after ex- 
plosion. Assuming s - 2 (i.e., B oc r'), it results in a magnetic 
field between 3.8 and 26.5 G at day 5 after explosion. 

If we use our approach, the spectral index, a - 1 .08 + 0.08 
(MarchiH et al. 2010; Brunthaler et al. 2010), is compatible with 
p ~ 3 or lower. However, using /3 - - 1 .43 + 0.05 (Marchili et al. 
120 lOi) and m ~ 0.89 (derived from a set of VLBI observations; 
Brunthaler et al. 1201 H and Brunthaler et al. in prep.), we find a 
self-consistent magnetic field of ~ 100 G for p -2.6 and s = 2.4 
(see Figs. |2] and lU, much larger than that reported in Marchili et 
al. (l20roi l. 

However, if p ~ 3, the magnetic field would be close to G, 
regardless of the value of s (in order to explain the spectral in- 
dex). Now, if we set 5 = 2, we obtain /3 1.3 for m = 0.89. 

This value is close to, but lower than, the observed one, and 
would increase if s would be slightly larger than 2. Indeed, the 
uncertainties in m, a, and /3 can still make possible s - 2 for 
p ~ 3.ln any case, the magnetic field in the emitting region can 
be arbitrarily small if p ~ 3, and we cannot favor neither this 
possibility nor the estimate of Bq ~ 100 G obtained for p - 2.6. 

In table l5.1l (row 5) we summarize the values of s, p, and Bq 
discussed for this supernova. 

5.5. Other RSNe 

In the cases of RSNe where only the radio light curves are avail- 
able, it is still possible to infer some information on magnetic 
fields and density structure of the CSM and/or ejecta, although 
with several additional assumptions. In this section, we study 
two cases which we consider interesting compared to other more 
typical RSNe. 

5.5.1 . Radio transient in M 82 

The discovery of a new transient in M 82 has been recently re- 
ported in Muxlow et al. (120091) . and a light curve with a prac- 



tically constant flux density has been reported in Muxlow et al. 
(120101) . with an spectral index of ~ 0.7. Indeed, looking at their 
Fig. 2, the flux density at 1.6 GHz seems to be slightly increas- 
ing. If this transient in the starburst galaxy M 82 is a supernova, 
it would be a so special case, since yS ~ 0. It is not possible to 
obtain such value of /3, unless s <2 (see Eq.|5]l, since the highest 
value of m is 1 and p is assumed to be larger than 1 . Indeed, from 
Eq.|5]we obtain s ~ 0.6, for p = 2, and s ~ 1.3, for the extreme 
case p - I. Therefore, a plain light curve is a strong evidence of 
a CSM density profile much shallower than the canonical case 
5 = 2. In any case, another condition for /5 ~ 0, regardless of 
the strength of the magnetic field, is that m ~ 1 (see Fig. [TJ. 
Therefore, two clear conclusions can be extracted for this tran- 
sient, provided it is a supernova: 1) the index of the CSM density 
profile is s < 2 and 2) the deceleration index must be m ~ I. 
Both conclusions imply that the index of the ejecta density pro- 
file, n, must be very large (n = 20, or even higher, see Eq.[T]). In 
regard of the spectral index, from Figs.[3]and|4]we conclude that 
the magnetic field would be up to Bq ~ 20 G, assuming s - 1.6 
and p -2, and lower for larger p. 

5.5.2. SN 2000ft 

Supernova SN 2000ft was discovered by Colina et al. (120011 ). 
Perez-Torres et al. (2009) presented an eight-year long radio 
monitoring of this supernova, located in the circumnuclear star- 
burst of NGC 7469 (a Luminous Infra-red galaxy, LIRG, at a dis- 
tance of 70 Mpc; Sanders et al. l2003l ). Perez- Torres et al. (|2009] ) 
followed the approach of Weiler et al. (120021 ) to fit the evolution 
of the radio light curves, using a standard value of i = 2 for 
the CSM. This analysis resulted in a value for the spectral index 
a - 1.27 and a power-law time decay index yS = -2.02. In addi- 
tion, they also needed to include a foreground absorber, likely an 
H II region, to account for the non-detection of radio emission 
at frequencies around and below 1 .7 GHz, in agreement with the 
observations reported by Alberdi et al. (l20Q6|). 

While the value of a reported for SN 2000ft is not surpris- 
ing, the value of /? is much larger (in absolute value) than those 
typically found in RSNe. From Eq.|5] it is possible to obtain val- 
ues of p similar to that of SN 2000ft if 5 > 2 (see Figs. [I]and|2]), 
although a low value of m (together with a large Bq) or a large 
value of p is also necessary to simultaneously explain the steep 
spectrum (see Fig. HJ. If i = 2, it is also possible to obtain a 
self-consistency between a and provided Z?o ~ 0, p = 3.54, 
and m - 0.75. 

In any case, we find that SN 2000ft should be a highly de- 
celerated supernova (m between 0.7 and 0.8), the CSM density 
index should be ^ = 2 or higher, and the energy distribution of 
the electrons must be quite steep {p -3 ox higher). 

6. Conclusions 

We have shown the impact of energy losses of relativistic elec- 
trons in RSNe, and how they affect the flux-density decay rate of 
the light curves in the optic ally- thin regime for different values 
of the magnetic fields and for different expansion curves. 

If the magnetic-field energy density and the acceleration ef- 
ficiency of the shock scale with the shock energy density, which 
is very likely the case for RSNe, we find that there is a tight 
relation between expansion index, m, spectral index, a, and (op- 
tically thin) flux-density decay index yS. 

This connection between expansion and flux-density evolu- 
tion in RSNe can be used to estimate the magnetic field of ob- 
served RSNe (Bo at a reference epoch) as well as its evolution 
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with time for an assumed CSM radial density profile and energy 
index, p, of the relativistic electrons. 

For a number of well observed RSNe (e.g., SN 1993J in 
M81), self-consistent solutions have been found for Bq, m, s, 
and p. A standard CSM density profile (i.e., s = 2) can ex- 
plain all observations, although evidences of non-standard val- 
ues of s ai-e found for SN 1986J and SN 1979C. The index 
of the relativistic electron population takes rather high values 
(p = 2.3 - 3.0) and the range of magnetic fields between all 
cases is large (Bq ~ 20- lOOG). These large magnetic fields im- 
ply effective amplification mechanisms in the radio-emitting re- 
gion, possibly related to plasma turbulence (see, e.g.. Gull 1 19731 
or Jun & Norman 1995 and references therein). 

Previous analyses of the radio light curves and expansion 
curves of these RSNe did not take into account the correct cou- 
pling between m, /3, and a for different magnetic fields. Some of 
the results previously reported for these supernovae could, there- 
fore, be internally inconsistent. 

The magnetic fields obtained with our approach are in sim- 
ilar to the equipartition magnetic fields. For SN 1979C and 
SN 1986 J, we obtain a range of self-consistent magnetic fields 
similar to those derived from equipartition with a lower acceler- 
ation efliciency for ions (i.e., low-to-intermediate values of k in 
Eq.[T3]). Additionally, for SN 1986J there is evidence of i < 2, 
provided the magnetic field is small. For SN 2008iz, either a very 
low magnetic field (with s ~ 2) or an extremely large magnetic 
field (with 5 > 2) are necessary to model the light curve, given 
the large flux-density decay rate (J3 = -1.43). For SN 1993J, we 
obtain a magnetic field similar to that reported in Fransson & 
Bjornsson (1998) and Marti- Vidal et al. (1201 Ibk although we 
use in our approach a subset of flux-density observations (and 
not the whole data set), to avoid possible biasing eff'ects coming 
from the ejecta opacity (Marti- Vidal et al. 201 IbV 

For the RSNe that will be detected in the future (the large 
sensitivity of the forthcoming radio observatories, like ALMA 
and SKA, will allow the detection and monitoring of many 
other RSNe), it will be necessary, in light of the results here 
reported, to study the connection between their expansion and 
flux-density evolution, in order to obtain self-consistent results 
for the CSM profile, the electron energy index, and the magnetic 
field, based on the observed spectral index, expansion curve, and 
flux-density decay index. 
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Appendix A: Derivation of tlie source function of relativistic electrons 

The source function is related to the acceleration of part of the electrons from the CSM, as they interact with the expanding supernova 
shock. If we assume that the acceleration efficiency scales with the energy density of the shock (see Fransson & Biornsson |1998l 
for a discussion of different possibilities and how they fit to the observations of SN 1993 J), the density of electrons instantaneously 
accelerated at a given time, t, will be 

ni-elOCHcsV^ (A.l) 

where V is the velocity of the expanding supernova shock and «cs is the density of the recently-shocked CSM, both quantities 
computed for the same time f. Since the radius of the shock is r oc f'", the number of relativistic electrons injected between r and 
r + dr is 

/ '■ 1 T/2 . _ „2 . 



Mel = F,,i n^, ^- j V'Anr' dr, (A.2) 

where s is the index of the CSM radial density profile, n^, and ro are the CSM density and shock radius at a given reference epoch, 
to, and Frei is the acceleration efficiency (or fraction of CSM electrons that are accelerated) at the same epoch Iq. 

These electrons are distributed according to oc E^^ (with E running from /imin = ni^ to infinity, being nie the electron mass). 
Therefore, the conservation of the number of electrons implies 

Mel = K I E-P dE. (A.3) 

Jm c- 

The factor K accounts for the normalization of the electron energy distribution. This factor is just 

K=^. (A.4) 
E P 

min 

Hence, the source function (i.e., the energy distribution of electrons shocked between r and r + dr) is 

S (E, r) = F,ei M,",L L E-P, (A.5) 

E P Vo 

mm 

which, in terms of time (given that r/ro - (t/to)"'), reduces to Eq.|7] 



Appendix B: Radiative and adiabatic energy loss vs. free-free loss in RSNe 

The rate of energy loss due to synchrotron radiation and adiabatic expansion, Ei and E^^, are given in Eq.[8]and|9] respectively. In 
regard of the energy loss due to free-free interactions with the CSM, we have (Pacholcz vk i 1 970i l 

£f~/l«es£, (B.l) 

where fi ~ (1 -8)x 10~'* in cgs units. The exact value depends on the level of ionization of the nuclei in the CSM (lower values of /i 
correspond to higher levels of ionization, which are expected in the shocked CSMQ Since an electron with energy E in interaction 
with a magnetic field B emits synchrotron radiation mostly at its critical frequency (given by v ~ ci B E^, where ci = 6.27 x 10'^ in 
cgs units, Pacholczvk l 19701 ), the ratio of radiative loss to free-free loss for electrons emitting at the critical frequency v is 



Ef f\ V ci no \to 



(B.2) 



where we have used the time evolution of B given in Eq. |2] On the other hand, the ratio of diabatic to free-free energy loss for 
electrons emiting at the same frequency is 

$ = (B.3) 

Ef Ji no r„ 

The ratios in Eqs. lB.2l and lB.3l evolve as power laws of time, whose indices depend on s and m. Therefore, for some combinations 
of s and m, the ratios will grow with time (and radiative and adiabatic losses will dominate over free-free losses), but for other 
combinations the ratios will decrease with time, and free-free losses may be comparable to the other contributions at late times. 



* There is an additional contribution to Ef that can be written as a modifying factor of /[ , which depends on log E. We have neglected this small 
correction. 
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In the case of radiative vs. free-free losses, the time index in Eq. lB.2l is positive if 
6 

m > , (B.4) 

6 + s 

which implies values of m larger than 0.75 for the canonical case s -2, although slightly higher values of m for lower s (for instance, 
m > 0.78 if s - 1.6). However, even if m is lower than these values (so the ratio Ei-/Ef decreases with time), radiative losses will 
still be higher than free-free losses for the times and observing frequencies of interest. For instance, with an initial magnetic field 
as low as Bo = lOG (at day 5 after explosion!) and a CSM density as large as no = 10^ cm"-', the ratio is ~ 15 at day 1000 after 
explosion, for m = 0.70, observing at 5 GHz. 

In the case of the ratio of adiabatic losses to free-free losses, it increases with time if m > 1/s. This relation allows expansion 
indices as low as 0.63 if s = 1.6, and even lower values for higher s. Even in the (so special) cases where we would find m < l/s, 
adiabatic losses would still be higher than free-free losses for the times and observing frequencies of interest. For instance, if 
m = 0.59 (i.e., the minimum possible value compatible with the Chevalier model for s = 1.6) and no is as large as 10^ cm"^, this 
ratio is ~ 10 at day 1000 after explosion observing at 5 GHz. 

In any case, for strongly decelerated RSNe (i.e., with low m) free-free losses might not be completely negligible, depending also 
on the CSM density (higher density implies larger free-free energy losses). In these special cases, the connection between /3 and m 
may not only depend on s, p, and Bq, but also on no. Hence, Eq.|6]will have to be individually integrated for each of these cases. 



Appendix C: Integral solution of the continuity equation in energy space 

The energy loss of an electron with energy E at time t is given in Eq. [TO] Here we rewrite the equation, 

2 . I t Y . E 

where p - 2{m - s) - 2 and Bref is the magnetic field at a reference epoch fref (we do not use the subindex 0, to distinguish the 
reference epoch from the initial time of integration, fo, see below). The solution to this equation is 

(1 +p-m)t"'Eo 

E{t) = 2 , (C.2) 

(1 + p _ ^)fm + 2/3 C2Bl^EotZit^^Pt^ - tl^^n 

where we have computed the integration constant by assuming E - Eo s.t t - Iq (i.e., the initial time of integration). We now write 
£0 in terms of E: 

(m - 1 - p)E 



Eo = ^ " , . —■ (C.3) 



(m - 1 - p)(to/tr + 2/3 C2Bl.Et-Wnto/tr - C") 



At a given time, fo, the electrons being shocked (let us call them N'^iEo)) are distributed as S (fo)£'g'' (where S (fo) is given in Eq. 
|7]|. Since the number of electrons is conserved, at a later time, f, these electrons (i.e., not all the electrons, but just those shocked at 
time fo) will have the distribution N'^{E) given by 

N°(E)dE = N''(Eo)dEo ^ N^E) = n\Eo)^. (C.4) 

dE 



Applying Eq. lC.3l and its derivative, to Eq. lC.4l we obtain 

N^iE) . Sir)im-l-pf-^E-^ir/rr ^^ ^^ 

((m - 1 - p)(T/f)'« + 2/3 C2Bl^ErP{{T/ty" tP^^ - t'+p)) " 

where fo has been replaced by t, which can take any value in the evolution time of the supernova. A similar approach was described 
in Pacholczyk ( |197()1 ) (see his Sect. 6.3), although a constant and homogeneous magnetic field was used. We notice that Eq. lC.5l is 
physically meaningful only when the power in the denominator is that of a positive number. Let us call tf the value of t such that 
the denominator of Eq. |C3] vanishes. For a given supernova, this value depends on E and f. For larger f, all the electrons shocked 
at time r will have energies below E. In other words, no electrons shocked at time tf (and earlier times) contribute to the electron 
distribution at time f for energies above E. 

It is now straightwforward to conclude that the total number of electrons at time f and energy E will just be the addition of all 
the (evolving) contributions of the source function between the beginning of the expansion, fo, and f (these contributions are given 
by Eq. lC.5l ). The resulting integral is 

NiE, f) . r Sir)im-l-pf-PE-Pirltrdr 

-J'F {(m-\- p)(T/f)"' + 2/3 C2B]^^Et^P^{(TltT tP^^ - t^^p)) 
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Appendix D: Why are the light curves flatter if we consider radiative cooling? 

From Eq. |3] we can approximate the value of yS by assuming that, at any time, the electron energy distribution does not differ so 
much from the canonical one, N oc . Under such an assumption (valid for a rough qualitative discussion), it can be shown that 

^ d\n{t) I 2 B N 

The term with BjB, which is negative and is related to the different electron emissivities under different magnetic fields, has the 
same effects on (3 either if radiative cooling is considered or not. In regard of the term N IN, which is positive and is related to the 
different total emissivity for different number of electrons, when radiative cooling is not considered it takes the form 



where c{E, t) - -Ve{NE) and E is given in Eq.fTO] The function c(E, f), which indirectly depends on A^, approaches as time goes 
by. Therefore, the numerator, N{E, t), if cooling is considered, will tend to that of the cooling-free case, while the denominator, 
N{E, t), will always be smaller in the case with cooling. Hence, 



From this equation, it is easy to see that the positive contribution to provided by N/N (see Eq. lD.lb is larger if radiative cooling 
is considered. The corresponding value of /3 is, therefore, closer to (the light curve is flatter). 

The inequaUty in Eq. lD.2l mav have important observational effects until very late times (several years), when we will have 



N(E,t) S(E,t) 



NiEj) ['s(E,t')d/ 



where S (E, t) is the source function given in Eq.[7] However, if radiative cooling is considered, we have instead 





(D.2) 




Even then, the average value of p in the whole optically-thin part of the light curve may still differ from that of the cooling-free 
case. 



